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Introduction 

The TSPT (Time Series Product Tool) software was custom-designed for NASA to rapidly 
create and display single-band and band-combination time series, such as NDVI (Normalized 
Difference Vegetation Index) images, for wide-area crop surveillance and other time-critical 
applications. The TSPT, developed in MATLAB®, allows users to create and display various 
MODIS (Moderate Resolution Imaging Spectroradiometer) or simulated VIIRS 
(Visible/Infrared Imager/Radiometer Suite) products as single images, as time series plots at a 
selected location, or as temporally processed image videos. 

MODIS is ideal for monitoring large crop areas because of its high temporal revisit rate (twice 
daily), its near global coverage, and its relatively small ground sample distance (250 m). 
Because MODIS imagery is acquired daily, rapid changes in vegetative health can potentially 
be detected. With the TSPT, noise removal and temporal processing techniques allow users to 
create low-noise, daily time series plots and image videos. MODIS metadata is used to find 
and optionally to remove bad, cloudy, and suspect data. The TSPT features the important 
capability of fusing data from the MODIS instruments onboard the Aqua and Terra satellites, 
which nearly doubles the effective temporal resolution. 

The TSPT GUI (graphical user interface) provides an interactive environment for crafting 
“what-if ' 5 scenarios by enabling a user to repeat product generation using different settings and 
thresholds. The TSPT Application Programming Interface provides more fine-tuned control of 
product generation, allowing experienced programmers to bypass the GUI to create more user- 
specific output products, such as comparison time plots or images. 

This type of time series analysis tool for remotely sensed imagery could be the basis of a 
large-area vegetation surveillance system. The TSPT has been used to generate NDVI time 
series to monitor crop phenology in California and Argentina and to monitor forest health in 
an area of southeast Mississippi following Hurricane Katrina. 


TSPT Data Flow 


Inputs 

• Science Datasets 

• Sensor Zenith Angle 

• Cloud & Quality Data 


Inputs... 

may be obtained 
from DAACs or other 
sources 

may be subsetted 
and/or reprojected 
with tools such as 
the MODIS 
Reprojection Tool 
(http://edcdaac. usgs. 
gov/landdaac/tools/ 
modis/index. asp) 
may be simulated 
with Application 
Research Toolbox 
(Ross et a!., 2006) 



\ Visualization Products 

• 1-D time series plots given a latitude, longitude, and 


date range 

• 2-D images given a region and a date 
^ 3-D time series videos given a region and a date range j 


Typical MODIS Input Datasets 


MOD02 Planetary Reflectance (Swath) 
MOD35 48-bit Cloud Mask 
MOD03 Sensor Zenith Angle 


MOD09 Surface Reflectance (Tile) 

MOD09 Metadata including Cloud 
Information 

MODMGGAD Sensor Zenith Angle 



Fused Aqua 
and Terra 


For each pixel, use 0 V , cloud 
data, and other metadata to 
choose between Aqua and Terra 

• PER PIXEL: If Aqua is 
chosen, use Aqua NDVI, 0 V , 
and cloud data 

• PER PIXEL: If Terra is 
chosen, use Terra NDVI, 0 V , 
and cloud data 
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Ideal Fused 
Aqua &Terra 
NDVI 


Fused 0 V 
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Cloud Mask 






Temporally Processed 
NDVI 

• Input can be from Aqua, Terra, or Fused 
Aqua & Terra NDVI. Similarly, for 
simulated VIIRS, the input could be VIIRS 
1030, VIIRS 1330, or Fused 1030 & 1330 
VIIRS. 

• Processing is performed in time dimension 
only. 

• Although handled internally as 3-D arrays, 
the processed outputs are 2-D raster grids 
(one for each day). 



Sheely Farm (2003) 



2003-08-14, Aqua & Terra 

Cloudy: 3.7%, Mixed: 0.4%, < Ideal: 1.9%. Mean SZA: 16.8 
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Color-coded Sheely 
Farm Crop Map 
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Fused and temporally processed NDVI from MOD02 
(Planetary Reflectance), MCD03 (includes 6 V ), and MCD35 
(Cloud Mask). Planted cotton and pima fields stand out from 
other Sheely Farm crop fields (8-14-2003). 


Fused Daily NDVI vs. MODI 3 16-Day NDVI 


Cotton Field 5-3 Aqua & Terra 
Latitude 36.215, Longitude -119.926 
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Fused and temporally 
processed NDVI using 
MCD09 and MCDMGGAD, 
with MCD13 16-day Aqua 
and Terra NDVIs shown for 
comparison. 

Note that the daily NDVI 
shows a drop off more than 
10 days before MODI 3. 
This difference would be 
critical in early detection of 
sudden threats to crop 
vi health . 


Katrina (2005) 


fbn right, fused and temporally processed Aqua and 
Terra NDVI from MCD02, MCD03, and MCD35 before 
the storm on 8-24-2005 and after the storm on 9-21- 
2005. The eye passed over Stennis Space Center, 
outlined in black, adjacent to the Mississippi-Louisiana 
\border. Bay Saint Louis is on right. 


NDVI plot at 
latitude 30.3458 
and longitude 
-89.6735, a few 
miles southwest 
of Stennis. 


NDVI Aqua-Terra 

Latitude 30.3458. Longitude -89.6735 




Before 


3M5-0&-74 Aqua & T«nri» 

Cloudy $5.3%. MUed: 12.4%. < Ideal 0.9*. Mean SZA 29.1 




4 


v 



After 


Summary 

• The TSPT provides a wide range of multitemporal analysis and visualization tools for 
detecting rapid changes in vegetation vigor via remote sensing imagery. 

• The TSPT makes judicious use of available metadata to locate and optionally to eliminate 
bad, cloudy, or otherwise suspect pixels. 

• The TSPT enables fusion of cross-calibrated sensors. 

Future Development 

• To support rapid prototyping efforts, integration with the Application Research 
Toolbox, used to simulate future data sources (Ross et al., 2006), is underway. 

• To provide easier integration with other tools, portions of the TSPT are being 
extracted as separate, stand-alone modules able to input and output common data file 
formats, such as GeoTIFF, ENVI® BSQ (band sequential), and HDF (Hierarchical 
Data Format). 
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